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Abstract 

In many calculations of the two-photon decay of hadronic molecules, the decay matrix 
element is estimated using the wave function at the origin prescription, in analogy to the 
two-photon decay of parapositronium. We question the applicability of this procedure to the 
two-photon decay of hadronic molecules for it introduces an uncontrolled model dependence 
into the calculation. As an alternative approach, we propose an explicit evaluation of the 
hadron loop. For shallow bound states, this can be done as an expansion in powers of the 
range of the molecule binding force 1//3. In the leading order one gets the well-known point- 
like limit answer. We estimate, in a self-consistent and gauge invariant way, the leading range 
corrections for the two-photon decay width of weakly bound hadronic molecules emerging 
from kaon loops. We find them to be small, of order C(me//3 2 ), where m and e denote the 
mass of the constituents and the binding energy, respectively. The role of possible short- 
ranged operators and of the width of the scalars remains to be investigated. 
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1 Introduction 



Hadronic molecules are bound states of two hadrons held together by the strong interaction - 
clearly to be distinguished from the so-called hadronic atoms, where the two hadrons are bound 
by the Coulomb interaction. In the latter case the strong interaction only leads to a slight shift 
in the binding energies (and an additional width). Hadronic atoms can nowadays be produced in 
laboratories almost routinely. Hadronic molecules, on the other hand, might well be part of the 
hadron spectrum but are not yet identified unambiguously. In recent years evidence has grown 
that a few of the large number of known scalar mesons might be of molecular character. For recent 
reviews on the meson spectrum, with emphasis on the heavy states, see Refs. pQEiE]. 

It was argued for many years that the studies of the two-photon decay of scalars could dis- 
tinguish among different scenarios for scalar meson structure. One of the most studied cases is 
that of the light scalar mesons a (980) and / (980) and indeed, the predictions of various models 
for these differ drastically. Assuming them to be qq states made of light quarks, one gets about 
1.3 to 1.8 keV for the /o(980) — > 77 width in the relativistic quark model [4j, while, under the 
ss assumption, the two-photon width of the /o(980) is calculated to be about 0.3 j- 0.5 keV [5]. 



Within the molecular model for scalars, the predictions vary from 0.2 keV in Ref. [6j*j to 0.6 keV 
in Ref. [9] and to 6 ke V in Ref. [TO] . In the present paper we demonstrate that the technique used 
in Refs. [9], [10] has a large theoretical uncertainty. We also show that a gauge-invariant treatment 
of the two-photon decay amplitude of the KK molecule yields the value of the 77 width for the 
/o(980) close to 0.2 keV. 

It is well known since long ago that the two-photon decay rate for the parapositronium is very 
well approximated by the product of the square of the wave function at the origin times the e + e~ 
annihilation rate at rest [TTl[T2"] . This was taken as a recipe by many authors and was applied also 
to calculate the two-photon decay rates of hadronic molecules [9j [10] . In this paper we argue that 
this procedure leads to wrong results. Instead we propose to calculate explicitly the hadron loops 
employing an expansion in the range of forces, 1/(3. Then the leading term assumes a point-like 
molecule vertex and the two-photon decay of a scalar meson is found to be 



77 



\ixJ \m s J 



2m\ 2 / ms 

arcsin 

ms J \2m 



ms = 2m — e, 



(1) 



where ms is the scalar meson mass, e denotes the binding energy of the hadronic molecule and 
m the mass of the constituents — for simplicity we only study systems with constituents of equal 
mass — and a = e 2 /47r denotes the fine-structure constant. The factor £ is different from 1, if 
not all constituents participate in the decay. For example, in case of the fo only the charged kaons 
couple to the photons in leading order and therefore ( 2 = 1/2. We also calculate the leading 
range corrections to the two-photon decay rate. They turned out to be suppressed by a factor 
me I ' (3 2 <C 1, since we focus on shallow bound states. 

The paper is organized as follows: in the next section we present some very general arguments, 
why the wave function at the origin cannot be used to calculate the decay of hadronic molecules. 
This will be demonstrated explicitly in the subsequent sections: in Section [3] we give the general 
formulae for the two-photon decay of bound states that allow us to investigate two limits: the 



^ tl In this paper the chiral unitary approach is used. That the scalars produced are to be interpreted as dynamically 
generated is shown in Ref. [7]. For a somewhat different view on this subject see Ref. |S]. 
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Figure 1: Illustration of the different parts relevant for the decay of a hadronic molecule into two 
photons. The vertex for the transition of the molecule into its constituents is denoted by r, the 
corresponding wave function by and the annihilation potential by A. Solid and wavy lines 
denote the propagation of the constituents of the molecule and of the photons, respectively. 



weak coupling limit — that leads to the wave function at the origin prescription — is discussed 
in chapter 13.11 and the limit of the point-like interactions in chapter 13.21 In Section H] the leading 
range corrections to the latter limit are calculated. We close with a summary and outlook. 

2 The relevant scales 

Before we go into details let us present some general arguments why the wave function at the 
origin should not be used to calculate the two-photon decay of hadronic molecules. The most 
obvious argument is that we simply do not know the wave function at the origin. In contrast to 
the parapositronium decay, the equations solved for hadronic molecules are not solved using the 
fundamental degrees of freedom. Instead one typically works with conveniently chosen interpo- 
lating fields — and this choice influences the short-range behavior of the molecule wave function. 
For the deuteron wave function this is to some extend discussed in Ref. [13]. Only the tail of 
the wave function is completely determined by the binding energy and is therefore known model 
independently. Our ignorance about the wave function at the origin translates into a large spread 
for predictions for the corresponding two-photon decay rate of, say, the light scalar mesons, from 
0.6 keV in Ref. [9] to 6 keV in ref. pffi 

The second argument is that any transition matrix element using the wave function at the 
origin meets certain problems with gauge invariance. In case of positronium this is a minor effect, 
since the violations are suppressed by at least one extra power in the fine structure constant a. 
In case of hadronic molecules, however, this violation might well be more severe. This will be 
discussed in some detail below. 

The third argument is that the hierarchy of scales in case of the decay of hadronic molecules 
is very different to that of positronium decay. The individual parts of the decay are illustrated 

# 2 Also the recent attempt to improve on the wave function at the origin formula presented in Ref. [14] is not 
a solution, for it suffers from the same ignorance and, in addition, leads to a violation of gauge invariance, as 
explained below. 
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in Fig. [TJ First of all there is the molecule vertex r for the decay of the molecule into its 



constituents — here two mesona# 3 l. Next come the two meson propagators. The final piece is 
the annihilation potential A, given by the photon-meson vertices and the intermediate meson 
propagator. Corresponding to the building blocks there are three scales relevant for the two- 
photon decay of the bound state. To begin with, there is the intrinsic scale rp of the vertex 
function r set by the dynamics of the bound-state formation. An additional scale r £ ~ 1/k 
appears due to presence of the bound state, where we defined the binding momentum k = ^/me, 
m is the mass of the molecular constituents. The third scale is given by the range of the annihilation 
potential. For a shallow bound state, the energy carried away by the individual photons is of the 
order of m. Consequently, the range of the annihilation is given by ta = 1/m. 

Let us consider parapositronium decay from the point of view of the hierarchy of the scales 
introduced above. We clearly deal with a nonrelativistic system with the binding energy e = 
a 2 m e /4, with m e denoting the electron mass. Note, it is the parameter k = y/m e e = am e /2 
that defines the long-range piece of the molecular wave function, which takes the form \l/(r) = 
a/ k 3 /tt exp (— kt). The vertex function depends only on the electron three-momentum p and is 
trivially related to the bound-state wave function [13] : 



r(p) = V2^(p 2 + ^(p), (2) 

with ms being the positronium mass. We denote wave functions in coordinate space by \I/ and 
their Fourier-transforms in momentum space by ip. An explicit calculation with the positronium 
wave function yields 

r(p) oc , (3) 

so that r r ~ 1/k in the positronium case. Finally, as discussed above, we have ta ~ l/m e . 
Therefore there is the hierarchy of scales 

Case A: ta -C r £ « rr- (4) 

Thus, in case of the decay of positronium, the annihilation process is well approximated as taking 
place at the origin and consequently the decay amplitude scales to an excellent approximation 
with the wave function at the origin. 

Quite an opposite situation takes place for molecular hadronic systems. Indeed, in this case 
the scale of the vertex function is defined by the range of binding forces 1/(3. If one deals with a 
loosely bound state formed by zero-radius forces (j3 — > oo) the hierarchy is 

Case B: r r < < r e . (5) 

Then annihilation process cannot be described with the wave function at the origin prescription. 

To see which case (case A or case B) is more adequate for hadronic molecules let us focus on 
the two-photon decay of the /o(980) as a kaon molecule. Then we have (3 ~ m p , where m p is the 
mass of the p-meson, the lightest meson participating in the meson exchange between kaons (there 
is no one-pion exchange between two pseudoscalars) , e < 0.1m, and, again, ~ 1/m. This leads 
to 

rr < r A < r e . (6) 



# 3 For simplicity we talk of mesons only for the constituents. Note that the reasoning does not need to be changed 
in the presence of fermions. 
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Figure 2: Graphical representation of Eq. (JTj) for the full meson-meson scattering amplitude. 



Comparing this to Eqs. (jlj) clearly shows that the decay of hadronic molecules calls for a very 
different treatment as compared to that for the decay of positronium. The question arises if it 
is at all possible to give a simple recipe to calculate such a two-photon decay of, say, the f . In 
this paper we argue that the corresponding decay amplitude is well approximated by a kaon-loop 
integral evaluated in the limit of a point-like decay vertex {(3 — > oo). We thus propose to work in 
the limit (j3j), as the zeroth approximation, and build finite-range corrections in powers of 1/(3 to 
the leading term. Naively one would expect such corrections of order of (m/(3) 2 which, in case of 
the f , turn out to be of the order of 40%. However, an explicit calculation presented below shows 
that the leading range corrections in 1/(3 only scale as (k//3) 2 which, in case of the f , turns out 
to be of the order of 1%. Therefore, in case of the fo the corrections to the point-like formula, 
Eq. ([T]), are expected to be at most of order (m//3) 4 ~ 15%. 



3 Bethe— Salpeter approach 

In this section we employ an explicitly gauge invariant approach based on the Bethe-Salpeter 
equation for the molecule vertex in order to illustrate in more detail the interplay of the various 
scales. For scattering amplitudes similar formalisms were discussed in Refs. [16j [TFJ HHj. The 
relevant equations in the two limiting cases A and B introduced above will appear as special cases 
of these general equations. 

Consider a Lorentz-covariant theory describing the meson-meson interaction via a potential 
V(p — k) which possesses the inverse interaction range (3. A priori no assumption needs to be made 
on the structure of V, however, to keep the expressions simple in this very general discussion we 
assume that there is no charge flow in the potential. As a consequence there will be no meson ex- 
change currents, when we include photons. This situation is naturally realized for potentials given 
by t-channel exchanges of neutral particles. This gives rise to the so-called ladder approximation 
for the scattering equation that we will refer to in the following for simplicity. 

In practice the just described restriction on the potential implies the omission of many di- 
agrams without a priori justification. However, the goal of this section is to demonstrate that 
in gauge invariant approaches self-energies get linked to scattering potentials. We will not draw 
any quantitative conclusions from the considerations in this section. In contrast to this, when we 
discuss the leading finite range corrections in Sec. 0] as well as Appendix [E] we do not need to 
make any additional assumptions and — to this order in the range of forces — the problem is 
solved exactly, next section the effect of charge exchange is explained in Appendix [El 
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Scattering of two mesons can be described by the equation (see Fig. [2]) 

T(p, q, P) = V(p-q)-iJ -0^S(k)S(k - P)V(p - k)T(k, q, P), (7) 

where P M is the total momentum of the bound state and p^ is the four-momentum of one of 
its constituents. The propagators given are solutions of the Dyson equation, presented in the 
graphical form in Fig. [3] — if we again assume V to refer to the emission and absorption of a 
neutral meson, we here work in the rainbow approximation 

S-\P) = S \p) - E(p), S (p) = - J-j, E(p) = -if ^S(k)V(p - k), (8) 

with mo being the bare meson mass. The physical meson mass m appears as the pole of the 
dressed propagator S(p). 

If there exists a bound-state with the mass m s , we may define the corresponding vertex function 
r(p, P) as the solution of a homogeneous Bethe-Salpeter equation 

r(p, p) = -ij j0pS{k)S{k - P)V( P - k)r{k, p), (9) 

which is to be evaluated at P 2 = m 2 s . The bound-state vertex is normalized through the condition 

~ 1 1 (0 r2(fc ' P ^w ms{k ~ p) = 2P ^ (10) 

which relates the vertex r(p, P) to the bound-state mass. 

To describe radiative processes one should first define the dressed photon emission vertex for 
a meson. In the absence of charge flow in the potential V this is (see Fig. HJ) 

/d 4 k 
^-^T(p, k, q)S(k)S(k - q)v^\k, q) 

= vf{p,q)-i J -0^V(p-k)S(k)S(k- q )v,(k,q), 

where 

v^(p,q) = (2p-q) ll (12) 

and g M and p^ are the emitted photon and the emitting meson momenta, respectively. As follows 
from Eqs. (JTj) and ffTTT) . the dressed vertex v^(p,q) obeys the Ward identity, 

q^{p,q)=S-\p)-S-\p-q). (13) 

The two-photon decay amplitude for the bound state can now be evaluated with the help of 
the diagrams depicted in Fig. and with the dressed vertices and propagators involved (notice 
that the seagull vertex in Fig. [5(c) need not be dressed since, due to the Bethe-Salpeter Eq. 
the corresponding diagrams are already included into the definition of the scalar vertex). The 
resulting transition matrix element is 

W = eU* 2v W» v , (14) 
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Figure 3: Graphical representation of the Dyson equations for the dressed meson propagator and 
for the meson self-energy, Eq. (jSJ). 




Figure 4: Graphical representation of Eq. (Hip for the dressed photon emission vertex. 



where 

d 4 k 
W) 



= e 2 j 7 — i r(k,P)v tM (k,q 1 )v u (k-q 1 ,q 2 )S(k-q 1 )S(k-P)S(k) 



+ e 2 J P>u(k, q 2 )v,(k - q 2 , qi )S(k - q 2 )S(k - P)S(k) (15) 

- 2^e 2 | -^r(k,P)S(k-P)S(k), 

with qi }2 and 2 being the four-momenta and the polarization vectors of the two photons. The 
quantity appears to be gauge invariant, 

W^qm = W» v q 2v = 0. (16) 

To show this one may use the Ward identity, Eq. ( TIB"]) , to write 

W^qf = e 2 J -0- i r(k,P)S(k-P)S(k-q 1 )v v (k-q 1 ,q 2 ) 

- e 2 J -0^r(k,P)S(k)S(k-P)v u (k- qi ,q 2 ) 

+ e 2 j -^r(k,P)S(k)S(k-P)v v (k,q 2 ) (17) 






k-P 



Figure 5: Diagrams contributing to the scalar decay amplitude. 



— e 



d A k 



r(k,P)S(k)S(k-q 2 )v u (k,q 2 ) 



(2k)*- 

2qi, I -0-r(k,P)S(k)S(k-P). 



Now, using the bound-state equation for r(k, P) and the second line of Eqs. (ITT]) one may write 
" ^ -r{k,P)S{k-P)S{k-q 1 )v v {k-q l) q 2 ) 

( 18 ) 



(2k) 



J ^y 4 r(l, P)S(l - P)S(l) (v„(l - q U q 2 ) - v^(l - Ql , q 2 )) 



The same manipulations applied to the fourth line of Eq. ( TTTj) lead to Eq. ( Tl6]) . 

We therefore see that it is necessary that the vertex function r and the photon-meson vertices 
are constructed consistently in order to get gauge invariant amplitudes. In other words, using in 
the expression for the decay amplitude the molecule wave function together with bare vertices and 
propagators, inevitably leads to the violation of gauge invariance. 

For the decay S — > 77 involving real photons, Eqs. ( Tl6l) imply that 

(19) 



iW = M(P 2 Mq% - g^( qi q 2 )}el^ P = q x + q 2 . 
Then, for the scalar of the mass 777,5, the total width of such a decay can be evaluated as 
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mi 



64vr 



|M(m 2 5 )p 



(20) 



where the identity of the photons in the final state is taken into account in the overall coefficient 
in Eq. ([201). 

Equation ( 115]) is still general (up to the absence of exchange currents) and we may study it in 
both limits: case A (Eq. PJ) as well as case B (Eq. ([5])). Note, in order to simplify Eq. (TI5]) we 
need to assume the coupling to be weak. In a situation, where case A holds for strong couplings, 
the full system of coupled equations needs to be solved. In the next subsections both limits are 
discussed individually. 



3.1 Case A in the weak coupling limit 

In case A, > r^. As we shall see, the decay width in this limit can be derived from the general 
expression of Eq. (|T5]) under the assumption of weak coupling. Then one may neglect the dressing 
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effects and self-energies altogether — in Eq. (1151) all propagators and vertices can be replaced by 
the bare ones. As outlined above, this necessarily implies a certain violation of gauge invariance, 
however, those violations are suppressed by at least one power in the coupling that is assumed to 
be small. In the limit considered, the typical momentum in the loop is very small and one may 
replace ko in the photon vertices as well as in the strong vertex f by m. Then one may write in 
the rest frame of the scalar (P M = (2m — e, 0)): 



/ ^h S (k)s(k-p) 

J 27T 



[2V 



where only the leading pole is kept and non-relativistic kinematics is used for the mesons with 
momentum k. A similar expression appears for q\ <-> q 2 . What remains to be evaluated now is 
the three-dimensional integral 



, ,2 f d 3 k /r J (2k - qi )^2k - P - q x ) 



W„ = -te\/-J j^gW) ^ ' + <" " »,« ~ ft) + 2^ ) , (22) 

where we used the relation ([2]) and defined k^ = (m,k). The term in parenthesis refers to the 
annihilation potential. By assumption we have Tp which translates into k <C q\ t 2- Under this 
condition one may neglect all k dependence in this term, which then reduces to the annihilation 
potential at rest, and pull it out of the integral. The remaining integral is nothing but the definition 
of the wave function at the origin (in coordinate space). These altogether yield gauge invariant 
answer ( Tl9i) for the amplitude, with 

M(P 2 ) = i(e^. (23) 

Thus one arrives at the following expression for the two-photon decay width for the limiting 

case A, 

m 2 

Note that the final answer is gauge invariant, but this is true only for the leading term in an 
expansion in the potential V for the transition matrix elements. As explained, if one wants 
to improve the accuracy of Eq. ([2^1) it is insufficient to just keep the momentum dependence 
in Eq. (122j) . as proposed in Ref. [H], but also the meson self-energies are to be kept explictly. 
Consequently, Eq. fl24l) should only be applied in the weak coupling limit. 

There exists a prescription to calculate the two-photon decay amplitude by contracting the 
on-shell decay amplitude with the bound-state wave function ip(k) (see, e.g. |20j): 

f d 3 k ^ r ^ ^ i 

Woe ——rip(k) W(K + (k)K~(-k) ^77) . (25) 



„2 



r 77 = 2C 2 ^l^(0)| 2 . (24) 



(27T)' 

Since gauge invariance is preserved for the on-shell amplitude W(K + (k)K~(— k) —>■ 77), then 
the full amplitude (125]) proves to be gauge invariant automatically. In the leading nonrelativistic 
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approximation the k dependence of W can be neglected, so that Eq. (|2"oT) is identical to Eq. (Ell) . 
However, in general Eq. (I25p violates energy conservation: in the cm. frame the use of the on- 

shell amplitude in Eq. (1251) implies that the kaon energies fcio = ^20 = 

+ m 2 , while energy 

conservation requires kio = k 2 o = rris/2. This problem is discussed in detail in Ref. [2~T] . 

Simple recipes to restore gauge invariance in the presence of non-trivial vertex functions 
through new contact diagrams with the derivatives of this vertex, successfully used for decays 
like (f> — > 7/0 [IS], fail, since the photons are not soft. As a result, gauge invariance, preserved 
in the point-like limit, appears broken already to order (m//3) 2 ~ 0.4 (see Appendix ICl for the 
details), where we used for illustration with m = rriK and (3 = m p the parameters relevant for the 
/ . In the previous section we showed that the inclusion of the scalar vertex structure in a gauge 
invariant way requires an accurate consideration of the dressed meson propagators and photon 
emission vertices. 

As stressed before, the approximations necessary to come to the wave function at the origin 
prescription in case of the positronium decay were justified, since only terms of higher orders in 
a need to be neglected (such corrections can be taken into account systematically, see |22j). In 
a strongly interacting system, where the couplings are typically of order unity or larger, these 
steps are not justified: they lead to uncontrolled results and potentially large violations of gauge 
invariance. 



3.2 Case B: The zero— radius interaction limit 

We now study the other limiting situation, case B (Eq. (jSJ)). In this limit we may assume the 
vertex function to be point like (J3 — > 00), which leads to a constant vertex function r(p, P) = gso 
for, say, the decay of the f into kaons. Then all dressing effects can be absorbed in coupling 
constants and masses and thus bare (in form!) vertices and propagators may be used (but for 
different reasons as compared to the previous subsection). 

Then the matrix element (fl9l) can be found from the set of diagrams depicted in Fig. [SI 



W a = (gsoe 2 



d 4 k e| ■ (2k - gi) e* ■ (2k - P - gQ 

(2tt) 4 ((k - qi ) 2 - m 2 )((k - P) 2 - m 2 )(k 2 - m 2 ) ' 

W b = W a (l <-> 2), (26) 
d 4 k 1 



W c = -2(g S0 e 2 (el 



(2vr) 4 ((jfe - P) 2 - m?)(k 2 - m 2 Y 



where m, as before, denotes the meson mass. 

The two-gamma decay of scalars can be viewed as a particular case of a more general situation 
of the S — > decays, studied in detail in Ref. [23], with the vector particle V also taken to be 
a photon. The details of the calculations are well-known, and can be found, e.g., in Refs. [2U ESI 
27J[28]. Notice that, although all integrals in Eq. (|26j) are logarithmically divergent, the sum 



W a + Wb + W c is finitd #4 l. Thus, adding these three yields for the amplitude M introduced in 



# 4 An elegant way of extracting the amplitude M, Eq. |27|) . by reading off a finite coefficient at a specific 
combination of the four-momenta in W a was suggested in Ref. 28 . The problem of convergence of the integrals 
(1261) was also studied in detail in Refs. l24l l29l [T51 



10 



Eq. (USD: 



M(P 2 



(27) 



27r 2 m 2 

with 1(b) being the loop integral function, 7(6) = I(a = 0, 6) (see, for example, Refs. [24"l 128] for 
the definition of I (a, 6)), where b = mj^/m 2 , 



1(b) 



dz 



l-z 



dy 



yz 



1 — byz 



(28) 



The analytic expression for 1(b) takes the form: 



1(b) 



1 2 
26 + P 

1 1 

"26 ~ 262 



Vb 

arcsm — 
2 



In 



Vb + Vb^I 



1 71 



6 < 4 



6 > 4. 



(29) 



^6^4 

Finally, using Eqs. (120]) and fl27|) together, one arrives at the decay width 
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c 



167r 3 m 4 



\ v m z 



(30) 



with the only unknown parameter being the coupling constant g$o- For a loosely bound system 
with P 2 = (2m — e) 2 ~ 4m(m — e), e m, the condition ffTUl) gives the relation between the 
coupling constant gso and the molecule binding energy e [15J, 



9so 

47T 



32mvme. 



(31) 



Inserting Eq. (J3T]) into Eq. (J3QJ gives Eq. 

For ms = 980 MeV and m = 495 MeV, which translates to £ = 10 MeV, one arrives at the 
prediction 

T 77 = 0.22 keV, (32) 

for the two-photon decay of the scalar /o(980), which we refer to as the point-like model prediction. 
In the following we shall derive an estimate for the accuracy of this result. 



4 Leading range corrections 

In the previous subsection we investigated the limiting case of [3 — > oo. In this chapter we 
derive the leading corrections that emerge from finite values of (3 — we shall calculate the leading 
corrections in 1/(3. This should provide a valuable insight into how accurate the formulae of the 
previous chapter should be expected to be. 

For this we use a simple covariant model complying with the requirements of the previous 
chapter and thus providing a gauge invariant description of the two-photon radiative decay of a 
non-point- like molecular state. 
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We start from an effective meson interaction Lagrangian which is responsible for the point-like 
scalar formation and supply it with an extra momentum-dependent self-interaction: 



L** = + ^ [Wfa)]* ■ (33) 

The form of the Lagrangian (1331) is chosen such that, after inclusion of the e.m. field, it does not 
give rise to extra meson-photon vertices. Indeed, since the Lagrangian (1331) is written completely 
in terms of the real field (p^<p, the standard substitution <9 M — > — ieA^ does not touch it. As 
a result, the set of diagrams contributing to the molecule decay to two photons is not modified 
and is still exhausted with the three diagrams depicted in Fig. [5j Additional terms that arise 
from possible charge exchanges just lead to more complicated expressions but do not alter the 
conclusions. This is discussed in detail in Appendix [Ej The theory described by the Lagrangian 
fl33|) can be renormalized to the given order 1//3 2 . We present the necessary details in Appendix 
IA1 and briefly summarize the results here. 

The effective meson-meson interaction given rise by the Lagrangian (I3"3"j) is 

V(p-k) = \ 1 + ^(p-k) 2 . (34) 

Note that, in addition to the two terms given in Eq. (13~4|) also a term that scales as A 2 (s//3 2 ) 
emerges from Eq. (1531 . where yfs = E cm . However, since we shall work at the fixed s = mf , this 
term can be absorbed into Ai. The dressed meson propagator and the dressed photon emission 
vertex are 

S(p) = — 2 »> ^(P>9) = Z ^( 2 P- Q)it + v»(p,q), (35) 

p z — m 

where the renormalization factor Z and the explicit expression for v^p, q) are given in Appendix 
IA1 We have 5 M (p, q)q fl = 0, so that the Ward identity (TT3T) is preserved. From now onwards we 
stick to the renormalized, physical, value of the mass m. Besides that, v^p, q) does not contribute 
to the radiative 77 decay under consideration since (ve*) q 2 =0 = 0, with e* being the photon 
polarization vector. 

We turn now to the Bethe-Salpeter Eq. for a loosely bound system. One can check that, 
to order m 2 / f3 2 and ^Je/m, the Bethe-Salpeter Eq. ((9]) is satisfied with the vertex function (see 
Appendix lAl for details): 

-1 f, , Mp-P)\ t A 2 



r(p,P) = Z-'g s {l+ C (32 J j, (36) 
and the normalization condition ffTUl) gives (see Appendix [B] for the details): 

7§ = 32myMe (l + , (37) 



4tt v V P 



which, as (3 — > 00, reproduces the relation fl3T|) obtained in the limit of the zero-range interaction. 

In the weak coupling limit that we focus on here, the bound state formation should be controlled 
by non-relativistic momenta. As a consequence g e g, the effective coupling constant of the bound 
state to its constituents, should have corrections at most of the order of me/ (3 2 [3Ql [31], for the 
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scale m 2 does not appear in non-relativistic equations. To recover this result we need to use 
Eq. (pjj) at the bound-state pole, P = P with Pq = m 2 , and for on-mass-shell mesons, p = p 
with p\ = (p — Pq) 2 = m 2 , to get 



(3* 



where the factor Z 2 was put according to the rules of the LSZ reduction formula. The scaling of 
the corrections in Eq. f l38l) is in line with the estimates of Refs. [301 EH ESI E2] • Here we used that, 
for the given kinematics, p(p — P) = (2m 2 — m 2 )/2. 

The general form of the matrix element is given in Eq. (Tl9|) . Following the method proposed in 
Ref. [28] (see Appendix [D] for an alternative method) we notice that only the diagrams (a) and (b) 
in Fig. [5] give rise to the structure q\ v q2^ in the transition matrix element (fl9|) . Moreover, these 
two diagrams give the same contribution to W, so it is sufficient to consider only one of them, 

W; v = e 2 (g s j 4^r(k, P)S(k - P)v v {k - q u q 2 )S(k - qi )v,(k, Ql )S(k) , (39) 

and to read off the coefficient at the structure q\ v q2^ which appears after the introduction of the 
Feynman parameters and shifting the integration variable [28]. Notice that in this structure the 
Z-factors coming from propagators, from e.m. vertices, and from the norm of the scalar vertex 
cancel against each other, so that 

W a =c 2 f d 4 k (2k- qi ),{2k-P- qi ) u {l + ik{k-P)/p 2 ) 
^ ^ 9S J (2tt) 4 ((k- qi ) 2 -m 2 )((k- P) 2 -m 2 )(k 2 -m 2 ) ' { } 

The corresponding loop integral is finite and the result reads: 

M(m|) = {ml) + £^M« (m|) , (41) 
where M^ ' is given by the point-like result, Eq. (12 7p . whereas M^ 1 ' takes the form: 

M«K) = -C-^I», 6 = 4 > ( 42 ) 
2n z m z m A 



with 



/vl • r . . . . 1 - y - z + 2yz 

(y + z — 8yz) m(l — yzb) + 



o 



1 — yzb 

Thus, up to order 1//3 2 , the integrals 1(b) and 1(b) enter the two-photon decay amplitude in 
the combination (1(b) + 1(b)) /1(b): 

•(=1^(1 + ^) , (44) 

where the factor 1 + £(m 2 / (3 2 ) comes from gs (see Eq. ( 1371) ) and 1(b) +£(m 2 / f3 2 )I(b) appears from 
the decay diagrams. The integral 1(b) can be calculated analytically. The result reads: 

1(b) = (l - |) Kb) . (45) 
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Figure 6: Dependence of the width T 77 , defined in Eq. ([T]), on the value of the binding energy. 
So, for the total width, we arrive at an extremely simple formula: 

r 77 (/3) = r 77 (i + e^^ 1 ) = r 77 (i + o (^) ) , (46) 

where r 77 is given by Eqs. (PQ) and fl32l) . and the corrections of order of m 2 /f5 2 cancel against each 
other. Contrary to Eq. ( |38l) . here the cancellation is unexpected and non-trivial: since the photons 
carry away an energy of the order of the mass m, their momenta are the same and therefore at 
least one of the particles in the meson loop has a typical momentum of the order of its mass. 
Consequently there is no justification for the use of non-relativistic kinematics in the evaluation 
of the two-photon decay of scalar mesons. 

Evaluation of the actual coefficient in front of the structure (me)/(3 2 would require making 
assumptions concerning the details of the molecule formation which are model-dependent, though, 
given a particular model of this type, it is straightforward to apply the technique of the present 
work to establish this coefficient which is expected to be of order unity (see also Appendix IBI). 

Equation fj46|) is the central result of this work, for it shows that the predictions derived for the 
limit a point-like interaction should be quite accurate. If one assumes the coefficient £ also to take 
its natural value of order unity (Ai ~ A 2 ), we find that the leading range corrections to Eq. f[52"j) 
should be of the order of me//3 2 , which translates into a few percent in the decay amplitude. 
Therefore the accuracy of Eq. ([T|) should be given by the sub-leading range corrections that are 
expected to be of the order of (m//3) 4 , which is about 15% for the case of the fo- 

Another source of uncertainty for any prediction based on Eq. ([1]) is our ignorance on the true 
binding energy. To investigate this point, in Fig. El the dependence of Eq. (01) on e is shown, again 
using for illustration the parameters relevant for the /o, namely m = m^. As one can see, the 
dependence on e is quite moderate, once the binding energy exceeds 5 MeV. Therefore, even if e 
is varied between 5 and 20 MeV around 10 MeV — the typical value used above — the predicted 
two-photon width changes by less than 0.05 keV. 
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However, one should be aware of the following important disclaimers: 

• most of the hadrons — including the f — are unstable. Thus the concept of vertex function 
and binding energy is not well defined for those, and one should employ a multi-channel 
Bethe-Salpeter formalism. The quantity that should replace the bound-state vertex r(p, P) 
in all the formulae given above is the multi-channel t-matrix. The proof of gauge invariance 
proceeds along the lines similar to those given in Section [3j In the molecular case, for the 
energies around the KK threshold (and far away from the inelastic thresholds) the amplitude 
in the KK channel can be written in the scattering length approximation with the complex 
KK scattering length: 

a KR = - , , k 2 > 0. (47) 

In the limit k 2 — > 0, the coupling to inelastic channels is switched off, and, for K\ > 0, there 
is a bound state in the KK channel with e = K\/m. As shown in [23], the data on, say, ttti 
scattering near the KK threshold can be described in the scattering length approximation 
with k 2 around 50 -v- 100 MeV, and the ratio K\/k,2 of order unity. Thus, the hierarchy of 
scales in the case of unstable scalar is similar to the one considered above. The two-photon 
decay of an unstable scalar meson in the limit of point-like interactions was evaluated, for 
example, in Ref. [36] . More systematic studies of the problem of unstable particles will be 
subject of a future work. 

• Another issue is the possible presence of additional short-ranged operators, for example, of 
the type of vector meson exchanges studied in Ref. [6]. Estimates for these and their proper 
inclusion in the renormalization program also go beyond the scope of the present paper and 
will also be subject of a future work. 



Summary 

1. The ^(O) formula for slow particles annihilation does not work for the two-gamma decays of 
hadronic molecules. Not only are the results numerically uncontrolled, which is reflected in 
a wide spread of predictions for the decay f Q — > 77 width found in the literature, but there 
is also a potentially large violation of gauge invariance necessarily present in the derivation 
of the formula. 

2. Simple recipes to restore gauge invariance in the presence of non-trivial vertex functions 
through new contact diagrams with the derivatives of this vertex, successfully used for 
decays like — > 7/0 [IS], fail, since the photons are not soft. As a result, gauge invariance, 
preserved in the point-like limit, appears broken already to order 1//3 2 . We showed that 
the inclusion of the scalar vertex structure in a gauge invariant way requires an accurate 
consideration of the dressed meson propagators and photon emission vertices. 

3. For phenomenologically adequate values of e — 10 MeV and (3 ~ 0.8 GeV for the scalar 
meson / (980) our prediction for the two-photon width is 

= (0.22 ± 0.07) keV. (48) 



15 



Our result compares nicely with the experimental values for the 77 width of the light scalar 
/o(980) m 

r 77 (/ (980)) = 0.311°;?? keV, (49) 

and [35J 

r 77 (/ (980)) = 0.281°;°° keV. (50) 

The new experimental value 



r 77 (/ (980)) = 0.205l° : °83(stat)l°:^(syst) keV (51) 

gives an even better agreement. This clearly supports the molecular assignment for the 
/o(980). 



It has to be stressed that the uncertainty of our theoretical prediction (J48J) so far only includes 
our estimate of the possible influence of the structure of the vertex function for the scalar meson 
(about 15% for the amplitude). Neither was the possible influence of the finite width included 
nor possible additional terms from shorter ranged transitions. Both will be subject of future 
investigations. 

It should be emphasized that the main goal of our study was to quantify the effect of range 
corrections to the two-photon decay of hadronic molecules in a model independent way. Those 
we identified as parametrically suppressed compared to what is expected naively. This finding 
should not be changed by the inclusion of inelastic channels (like tttt in case of the /q). From 
this point of view our work is an additional justification for the use of, for example, the chiral 
unitary approach, for the calculation of the two-photon decay of the light scalar mesons [6] . Here 
scalar mesons appear as hadronic molecules based on point-like interactions. On the other hand, 
in Ref. [6] the 7T7T channel is included in a coupled channel framework. It is also reassuring that 
the width calculated in this reference is consistent with our result, Eq. ( 1481 ). 
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A Renormalization of the model and the Bethe— Salpeter 
equation 

To renormalize the theory ( |33j) to the given order 1/ j3 2 we start from the interaction Lagrangian 

L**=\\itf<P? + ^^V)]' (52) 
We consider the meson mass operator then, 

S(p) = -x J -0^S(k)V(p-k), (53) 

and use the dimensional regularization scheme to make it finite. It is easy to see that £(p) can be 
written in the form 

S(p) = (1 - Z-^p 2 + 6m 2 , (54) 

where 

*=1 + ^(A + D. A=f- 1 > + h i£, (55) 

with D = 4 — e being the number of dimensions, fi and 7^ ~ 0.577 being an auxiliary mass 
parameter and the Euler constant, respectively. The physical meson mass is simply m 2 = Z(m,Q + 
5m 2 ), and the meson propagator takes the form given in Eq. fl35l) . It is also straightforward to 
evaluate, to the same order l/(3 2 and in the same regularization scheme, the photon emission 
vertex, 

Vfl (p, q) = (2p -I)*- if J0U V (P ~ k)S(k)S(k - q)(2k - q)„ (56) 
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to arrive at 
where 



Vp(p,q)=Z (2p-q) fl + v ll (p,q) 



(57) 



A 2 



12np 2 



A 1 4m 2 — q 2 ( Am 2 — q 



3 + 6 



q 1 



+ 



2\ 3/2 



q 



arctan 



Am 2 - q 2 



\q 2 Vv ~ (P9)?a»] 



(5f 

with the renormalization factor Z given in Eq. (15511 . This agrees with Eq. (135]) . 

Now, before we come to the Bethe-Salpeter equation, we introduce two auxiliary integrals, 



d 4 k 



1 



(2tt) 4 (k 2 - m 2 )((k - P) 2 - m 2 ) 
d 4 k k 2 



(59) 



(2tt) 4 (k 2 - m 2 )((k - P) 2 - m 2 ) ' 
which are divergent and, in the dimensional regularization scheme, take the form: 



h{P) 



(4tt 



1 2m 2 

A + I QR (P), I 2 (P) = ——\ + I 2R (P), 



(4tt) 



(60) 



where Iqr and I2R are finite. 

The Bethe-Salpeter equation is 



r(p,p) = -i 



^P(k, P)S(k)S(k - P) (\, + j 2 ( P -k) 



(61) 



where dressed kaon propagators should be used. In the leading order in l/(3 2 the scalar vertex 
can be found in the form: 

r(p,p) = gi + ^p(p-p), (62) 

with the coefficients g\ and g 2 satisfying the equations: 



AiA 2 ^i 



m 



01 = Xigilo(P) + 2^^/o(P) 77Z ^(A + 1) 



Xl92 



UP) 



X192P' 

2(3 2 



UP) 



^291 

P 2 



(An) 2 " ' (3 2 

92 = \29ih(P)- 

These yield the equation which defines the mass M of the bound state as (P 2 = M 2 
A 2 



HP), (63) 



(64) 



1 - A, /„(/') 



1m 2 1 
Al/o(P) W) 2{k + 1} + Al/ °( P ) /2 ( P ) - 2 P2XJ ^ P) + h{P) 



0. (65) 



We treat this equation perturbatively in 1//3 2 . Then, in the zeroth order, one has 



t- = Io(M ), 



(66) 
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where M is zero-order mass of the bound state. The divergent part of the integral Iq(M ) is 
absorbed into the coupling constant Ai. Then, in the next -to-leading order, 



1 2f 



771 1 

— (A + 1) + / 2 (P)- 7 P 2 / (P) 



(4vr 



4 



0, £ = A 2 /A 



i- 



(67) 



with the problem of renormalization solved similarly to Eq. f[66|) . As g 2 enters the vertex together 
with 1//3 2 , Eqs. fl64l) and (1661) together yield 



92 = £gi, 

and thus we arrive at the vertex function in the form 



r(p,p) = gi li + z 



ff 1 



(6* 



(69) 



which requires normalization. This is discussed in detail in Appendix | 



B Normalization of the vertex function 

Normalization of the vertex function is given by Eq. fflOl) . In the zeroth order in the 1 /ft 2 expansion 
this gives 

which, for a loosely bound state with e -C m, reproduces the relation f[3~Tl) with g S o = Q\- Let us 
go beyond the zeroth order now and include corrections oc l//? 2 . The form of r(p,P) given in 
Eq. (169]) obviously represents the first two terms in the successive expansion of the exact vertex 
function in the inverse powers of f3 2 . Such an expansion performed prior to taking integrals with 
r(j>, P) involved may explode if the rest of the integrand does not converge fast enough to suppress 
the contribution of the higher and higher powers of the loop momentum which appear in the 1//3 2 
expansion of r(p, P). This is obviously not the case for the loop integral (fTUl) and thus we face the 
problem of convergence of the normalization integral already to order 1//3 2 . Notice, however, that 
the solution of this problem is well-known — the normalization integral evaluated with the vertex 
function in the full form converges, and the expansion in l/(3 2 is to be performed afterwards. 
Building the full form of the vertex function would require resorting to a particular model of 
the molecule formation which we would like to avoid in our general consideration. Fortunately, 
for loosely bound states, the model-dependent contributions to the vertex function appear in 
higher orders in the e/m expansion. Indeed, if we substitute the vertex function (|69|) to the 
normalization condition (1101) and perform integration in d A k retaining only the terms of order 
\Je/m and neglecting all higher contributions, then the result can be expressed entirely through 
the integrals Iq(P) and Iz{P) defined in Appendix [A] computed to the same order, 

T 1 . 1 1 fF T , 2m 2 . 3m 2 m 2 [T ,„„ . 

(47r)^ Stt z 167r V m {4:7T) Z lbir z Ibir V m 

All divergent contributions disappear, as discussed in Appendix [A] and the result (|37|) is readily 
reproduced, where 

9s = Z 9l . (72) 
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Then, finally, the scalar vertex function takes the form of Eq. (1361) . 

As a cross-check of the results (I36p and (1371) we assume a particular form of the full scalar vertex 
compatible with the large-/? expansion fl36|) . As such we choose, for the sake of transparency, the 
form: 

r <.P,n = Z-> SS ff2 _^ p _ py (73) 

Then, introducing Feynman parameters and integrating out the four-momentum, one can rewrite 
the normalization condition (fTUl) in the form: 

9l 



2 f 1 zdz f(^ 2 /^ 2 ) x{l-z-jxQ 2 /^ 2 )) 



8tt 2 J Q 2 J [x 2 eQ 2 P 2 /{3 4 + x(l + £(zP 2 -m 2 )/[3 2 ) + l} 

where Q 2 = m 2 — z(l — z)P 2 . The integral in x can be easily evaluated and yields, to order 1/ (3 2 

f^ 2 ^ x(1-z-£xQ 2 /(2P 2 )) ^ 1 / zP 2 -m 2 

J X [x 2 eQ 2 P 2 /f3 A + x{l + i{zP 2 -m 2 )/(3 2 ) + 1] 3 ~ 2 l Z) \ K (3 2 

so that the relation ( 1711) reduces to 



(75) 



, P 2 -2m 2 ^ g 2 f 1 z(l-z) 

where the symmetry of the function z(l — z)/Q 2 with respect to the variable change z —>■ 1 — z was 
used. The remaining integral in z is specific for the point-like limit and it was evaluated before in 
order to derive the relation (l37l) . Therefore, one can rewrite ( |76l) in the form 

-(?))(«)„.,, - 

Thus the relation (|37|) is re-derived plus the first correction of the form (me)/j3 2 is established 
for the vertex function (1 73 p . As it was anticipated before, the term of order m 2 / (3 2 is model- 
independent and coincides with the one obtained in the simple approach described in the beginning 
of this appendix. 

The result of this appendix can be understood in the language of effective field theories. Indeed, 
all divergencies are to be absorbed into appropriate counter terms, however, there are no counter 
terms allowed that are non-analytic in e. Consequently all terms that scale as \fe are fixed model 
independently. 



C Failure of a simple recipe of gauge invariance restora- 
tion for a non— trivial vertex function F{k, P) 

In this appendix we demonstrate that gauge invariance appears broken to order m 2 / f3 2 for the 
naive attempts to consider a non-trivial vertex function r(k, P) without a self-consistent dressing 
of the photon emission vertices and the meson propagators. Thus we start from Eq. (fT5|) with 
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a non-trivial vertex function r(k, P) but with the bare photon emission vertices and meson 
propagators. This yields: 



W^ql = e 2 J 4^i[r(k + Qi, P) - r(k, P)]S(k)S(k - q 2 )(2k - q 2 ) v . (78) 

For r(k, P) ^const the difference in the square brackets does not vanish, and as a result also 
W^qi remains finite. A naive counting of powers of f3 shows that this difference scales as l/(3 2 
(see, for example, Eq. (1361) for the form of the vertex r(k, P)) and thus a simple trick with adding 
contact diagrams with the photon emission from the scalar vertex described by dPjdk^ could 
solve the problem to the given order 1//3 2 , so that gauge invariance breaking happens only in the 
order 1//3 4 which we neglect throughout the paper. In Ref. [15] this approach was successfully 
used for the decay 4> — > 7/0 — the interested reader can find the details and the discussion of the 
method, for example, in the aforementioned paper. In contrast to the 0-decay, here we have two 
identical photons in the final state. The requirement of symmetry of the amplitude with respect 
to the interchange of these leads to two single contact diagrams and one double contact diagram, 
the latter containing d 2 Pjdk^dk v in the scalar vertex. The contributions to be added to W^q^ 
take the form: 

d'k dr(k,p) 

S{k)S{k - q 2 ){2k - q 2 ) v , 











5W$ q ? = 


*7 



(2tt) 4 dfc" 
d 4 k dr{k,P) 



S(k)S(k-q 1 )(2k-q 1 ) fl , (79) 



(2tt) 4 dk v 
> 4 h Ffrih p\ 

S(k), (80) 



so that the resulting expression reads 

d 4 k 



W^ = e 2 J 



(27T) 



(2tt) 4 dk^dk 

dr(k,p) 



r(k + qi ,p)-r(k,p)- q l t- 



(81) 

x [(2k - q 2 ) u S(k - q 2 )S(k) - 2k„S 2 (k)] , 



where the integral coming from the double contact vertex was integrated by parts. In order 
to proceed we use the form of Eq. f!73|) for the vertex function — since we are only after the 
scaling behaviour of the remaining violation of gauge invariance, here we are free to work within 
a particular model. This choice ensures convergence of the loop integral and complies with the 
large-/? expansion ( l36i) . Then we get 

P(k + qi ,P)- P(k, P) - qf-^1 = Z-W ^ - ( ^ 2)]2 



dk» aa {(3 2 - f (fc + qi)(k- q 2 )] 2 {/3 2 - ik{k q 2 )] ' 

(82) 

which scales as 1//3 4 as (3 — > 00. However, in contrast to this the corresponding integral of Eq. ( |HTT) 



for W^q^ scales as 1/(3 . To see this observe that convergence to the integral is provided by the 
denominator with the consequence that k takes values of the order of j3 (for large (3). Therefore, 
the relevant estimate for Eq. (182]) is f3 2 (kqi) 2 / /3 6 ~ uj 2 //3 2 , where us is the typical energy of a 



22 



photon. We checked by an explicit calculation that this behaviour is valid for the integral AST]) as 
well. 

Therefore, even introducing the correction terms of Eqs. flHUl) . does not change the order where 
a violation of gauge invariance appears — it appears at order u 2 /f3 2 . In the given kinematics 
photons are not soft, u ~ m. Thus the simple prescription described in this appendix to cure the 
violation of gauge invariance does not improve the situation, since gauge invariance is still broken 
by the terms of the order of m 2 /f3 2 . 



D Alternative derivation of Eq. (1421) 



As a cross-check of gauge invariance, let us extract the amplitude ( ]42l) from the coefficient at 
the structure g^ v in Eq. f[T§l) . This is less trivial as the seagull diagram (Fig. 0(c)) contributes. 
This seagull has Z" 1 factor coming from the scalar vertex and Z 2 factor due to the two meson 
propagators. Because of this mismatch of Z-factors, in addition to the contribution giving the 
result (1421) . a divergent piece arises in W^, which comes from the leading term, of order (1//3 2 ) , 
in the scalar vertex r(p, P): 



- 2(tg^g s e 2 I (Z - 1) = -2(tg, u g s e 2 I ^^(A + 1), (83) 

where the expressions for the Z and A in the dimensional regularization scheme are given in 
Eq. fl55l) . To order l/(3 2 , one can make use of Eq. fl66|) and replace A 2 /o i n flHUD by 

X 2 I = ZXJq = C (84) 

In the meantime, another divergent contribution comes from the 1 / j3 2 term in r(k, P), which 
reads: 

2&g, u g s e 2 j- 2 —(A + l), (85) 

so that the two undesired divergent contributions to W^ v cancel against each other and the gauge- 
invariant formula (fl9l) is re-derived, with M(P 2 ) given by Eq. (f4~TT) . 



E Account for exchange currents 

The interaction Lagrangian (133|) does not give rise to extra kaon-photon vertices, in addition to 
those following from the kinetic part of the kaons Lagrangian. In this appendix we consider the 
possibility for these vertices to appear due to meson exchange currents. 

Let us introduce field doublet ip a , and have the interaction Lagrangian in the form 

Lint = ^A^+fV) 2 + ^- 2 [dMw)} 2 . (86) 
In momentum space, it gives rise to the four-point vertex of the form 

(ai + ^(p-^) 2 )^ (87) 
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where p and k are the momenta of the incoming and outgoing kaons, and a, (3 and p, A are the 
isospin indices of incoming and outgoing mesons, respectively. As before, the term linear in s was 
absorbed into Ai. 

Electromagnetic interaction is then found from the minimal substitution, 

p^^p^-eA,Q% (88) 

where Q — (1 + Ts)/2 is the charge operator. Thus two new kaon-photon vertices are generated 
in the order 1//3 2 : the contact single-photon vertex (q^ is the photon momentum), 



i^e 3kn r^ p (2k-2p + q)„ (89) 

and the double-photon vertex, 



2 ^(^-^3>3p)- (90) 

The dressed propagator is now 



S$ = S2S(p), s(p) = -^2i (91) 

with 

3X2fn 2 , . . . . 

z=l+ ww (h+1) - (92) 

The dressed photon emission vertex, 

7-1 



^M = Qp&P-Q)*-i^Q P aTp I ^S(k)S(k-q)(p-kf(2k-q) 
+ ^e 3kn r^J -0- A S(k)(2k-2p + q)„ 



(93) 



satisfies the Ward identity 

q^ = Q X p [S- 1 {p)-S-\p-q)l (94) 

with S(p) given by Eqs. (|9T|) and (I9~2l) . After simple algebraic transformations the photon emission 
vertex (1931) can be represented as 

^(p,g)=Q^- 1 (2p-g) M + ..., (95) 

where the ellipsis denotes the terms irrelevant to the 77 decay involving real photons. 
The Bethe-Salpeter equation is now 



I>, P) = -i J (04^^ P)TpS(k)S(k - P) (\, + J 2 (P~ A:) 2 ) • 



(96) 



The momentum dependence of the vertex and the normalization condition are the same as in 
Eqs. fl36|) and fl37j) . respectively, while the matrix structure of the vertex is either 

/>,P) oc-U" (97) 
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(d) 



(e) 



Figure 7: Extra diagrams contributing to the two-photon scalar decay amplitude due to exchange 
currents. 



in the isosinglet case, or 



r?(p,p) oc-^, 



1 

3' 



(98) 



in the isotriplet one. Consequently, in the former case, Eq. (1661) is replaced by 



whereas, in the latter case, it becomes 



T = 3/o(P), 



(99) 



(100) 



Thus one may have either a isosinglet or a isotriplet bound state, depending on the sign of Ax- 

Let us turn to the calculation of the two-gamma decay amplitude. The contributions to the 
decay amplitude proportional to q\ v q2^ are left intact by inclusion of the exchange currents and 
they are given by the graphs (a) and (b) of the Fig. [51 In the meantime, extra terms appear 
which contribute to the coefficient at the structure g^. Because of the mismatch of Z-factors, 
the divergent piece coming from the leading order term of -T"(p, P) appears, similarly to Eq. fl83l) . 
to be 

-e^^/oA^^A + i), (101) 

while the contribution from the 1//3 2 term in r^ip, P) remains the same as in the neutral exchange 
case — see Eq. ( 185]) . Thus, using the relations (1991) and fllOOp . one arrives at the conclusion that, 
in the isosinglet case, the graphs of Fig. [5] do not generate extra contributions proportional to g^ 
either while, in the isotriplet case, such an extra term reads: 



9 /- 2 

gse £ m 



(A + 1). 



(102) 



In addition to the graphs depicted in Fig. [5j there are extra contributions due to the presence 
of new contact vertices (1891 and (J9~0|) (see Fig. |7]). The single-photon vertex fl89l generates the 
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W $ = 9se 2 ^e 3kn r^r^ p Q^ 



contribution (see Fig. 0(d)) 

P 2 

x J ^y{2*y ms{k ~ p W q W q - q ^ 2q - 2k + qi ^ 2q ~ q ^ ( 103 ) 

+ (1 <-> 2, /U <-> z/), 

where the matrix structure of is given by either Eq. ( 1971) or by Eq. ( |98l) . The double-photon 
vertex ( 1901) gives rise to (see Fig. [7J(e)) 



Wg) = 2 Ms e 2 ^|rf (^ - t$t£) | 



dfc dg -5(A;)5(A;-P)5(g). (104) 



(2vr) 4 (2tt) 



Explicit calculations yield, for real photons, 

W$ + W$ = 0, (105) 

in the isosinglet case, and 



in the isotriplet case, so that this divergent term cancels against that given by Eq. (11021) . As a 
result, we conclude that Eq. ( 1461) holds also if exchange currents are included. 



26 



